B*B*p vertex from QCD sum rules 



Chun-Yu Cui, Yong-Lu Liu and Ming-Qiu Huang 

Department of Physics, National University of Defense Technology, Hunan 410073, China 

The form factors and the coupling constants in the B*B*p vertex are evaluated in 
the framework of three-point QCD sum rules. The correlation functions responsible 
for the form factors are evaluated considering contributions of both B* and p mesons 
as off-shell states. The obtained numerical results for the coupling constants are in 



agreement with light-cone QCD sum rules calculations. 
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I. INTRODUCTION 



In the last years, both experimental and theoretical studies on heavy mesons have 
received considerable attention. With growing datasets collected by the collaborations 
D ■ such as CDF, DO, CLEO, and the forthcoming SuperB and 



spectroscopy and decay of heavy flavor become more exciting [1 



HC, investigations of the 
4|. Therefore, the reliable 



determination of various characteristics, such as form factors and coupling constants are 
X* 1 needed. 

r*"* \ Suppression of heavy quarkonium has been considered for a long time as one of the 

^ \ most striking signatures for the quark-gluon plasma (QGP) formation j^J. It results from 
Tj- \ Debye screening of color force in the QGP. In this picture, due to larger abundances of 
^ ! color charges screening the interaction between the c and c quarks, the attractive cc po- 
tential responsible for the J/ip binding gets screened as the temperature of the medium 
increases in the QGP. (For a detailed review see Refs. 0, 0] and references therein.) Ex- 
tensive experimental efforts have been devoted to study this phenomenon at the Super 
Proton Synchrotron (SPS) at CERN and Relativistic Heavy Ion Collider (RHIC) at the 
Brookhaven National Laboratory. As the case of charmonium states in the QGP which are 
sensitive to the color screening effect, the study of bottomonium suppression in high energy 
heavy ion collisions can be used as a signature for the QGP as well [7]. In contrast to char- 
monium states, bottomonium states should be a cleaner probe of QGP due to its low cross 
section where the competing effects, which either reduce the yield Q or enhance it (91], are 
negligible. We expect the effects of the QGP on the absorption of T in ultra-relativistic 
heavy ion collisions at the LHC, which can provide an answer to this open question. On the 
other hand, other more conventional mechanisms based on J NAT) absorption by comoving 



hadrons have also been proposed as a possible explanation [10Ml2f| . In this way, one needs 
to understand the effects of T absorption in hadronic matter. It is known that tt and p are 
the dominant hadrons in ultra-relativistic heavy ion collisions. Since it is still difficult to 
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study strong interaction phenomena at non-perturbative regime using the QCD, the study 
of quarkonium absorption is generally performed in the framework of effective Lagrangian 
with meson exchange 0]. In the process, T and p produce the final states B* and B* by 
exchanging a B* meson. The calculation of T absorption cross sections needs information 
of the B*B* p interactions. To describe the strong interactions of the negative-parity heavy 
mesons with p meson, we employ an effective Lagrangian, which is constructed based on 
the chiral symmetry and described by the following 131 ] : 



C = WBBpTr 
+ igB*B* P Tr 



B^ d^B) P^ 



- 2f B * Bp e^Tr [B^ d p B* v - B*J d p B) d a P f 



B*J d u B*» ) P h 



4if B »B* P m B *Tr [(b$B*) {d»P v - d v P») 



where B and B* represent isospin doublets, P is the isospin triplet of the p meson. The 
B*B*p interactions are characterized by two indepentent coupling constants g B * B * p and 
Jb*b*p- It is necessary to know values of coupling constants with some precision. The 
choice of a lower or higher value may change the final cross section to some extent. 

Theoretically, the knowledge of the heavy-heavy-light mesons coupling constants in 
hadronic vertices are very important in estimating strength of hadron interactions when 
hadronic degrees of freedom are used. They are fundamental objects of low energy QCD. 
They may also play an important role in the formation of these possible molecular can- 
didates composed of two B mesons. However, such low-energy hadron interaction lie in 
a region which is very far away from the perturbative regime, precluding us to use the 
perturbative approach with the fundamental QCD Lagrangian. Therefore, we need some 
non-perturbative approaches, such as QCD sum rules(QCDSR) jl4]-ll6l|. to calculate the 
form factors. Besides, QCDSR at finite temperature illustrates mass shifts and width 
broadening 0, 0. Thus, it is expected that the form factor of B*B*p vertex may be 
sensitive to high temperatures and it is maybe another interesting task. 

In Ref. 0, the coupling constants g B * B * p and f B * B * p are estimated using QCD light 
cone sum rule (LCSR) method. In this article, the form factors and the coupling constants 
of the B*B*p vertex is calculated in the framework of the three-point QCDSR. We notice 
that in the case of D*D*p vertex, the form factor and coupling constant have been studied 
with three-point QCDSR using the effective Lagrangian based on SU(4) flavor symme- 
try {23)]. Different from above situations, we consider the effective Lagrangian, which is 
constructed based on the chiral symmetry. Herein, we use the same technique developed in 



the previous work for the evaluation of the couplings in the vetices D*Dn [19j,|20|], DDp [21] , 
D*D*tt 0, D*D*p 0, DDu 0, D*DK* (892) 0, D S DK* and B* sl B*K 0. 

This paper is organized as follows. In Sec. [Ill we give the details of QCDSR for the 
B*B*p vertex when both B* and p mesons are off-shell. Sec. IHIl is devoted to the numerical 
analysis and discussion. Additionally, the Appendix presents the formula of form factors. 
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II. THE SUM RULE FOR THE B*B*p VERTEX 

In this section, we give QCDSR for the form factors of the B*B*p vertex. The three- 
point function associated with the B*B* p vertex, for an off-shell B* meson, is given by 

rf;:W) = / d 4 xd 4 y ^'■ x e-^{Q\T{f;{x)jT\y)3Tm^), (2) 

where the interpolating currents are j p ( x ) = ^( x )lfj,d(x), j„*°(x) = d{x)^ v b{x) : and 
3a i x ) = u(x) , y a b(x) . The correlation function for an off-shell p meson is 

= J d'xd'y e*->e-**» (0|T{jf \x)jf (y)j** t( )}| > . (3) 

In the expressions, q = p'—p is the transferred momentum. There are fourteen independent 
Lorentz structures in the general expression for the vertices fl2]) and ([3]). We can write T^ ua 
in terms of the invariant amplitudes associated with each one of these tensor structures in 
the following form: 

T fva(PiP') = t i(p 2 >p' 2 , q 2 )g^p a + r 2 (p 2 ,p' 2 , q 2 )g m Pv + T 3 (p 2 ,p' 2 , q 2 )g va p P 

+T A (p 2 ,p' 2 , q 2 )g^p' a + T 5 (p 2 ,p' 2 } q 2 )g^ a p' u + T 6 (p 2 ,p' 2 , q 2 )g va p'^ 
+T 7 (p 2 ,p' 2 , q 2 )p^p v p a + T 8 (p 2 ,p' 2 , q 2 )p^p v p a + T g (p 2 ,p' 2 , q 2 )p^p' u p a 
+ri (p 2 ,p' 2 , q 2 )PpVvP a + T n (p 2 ,p' 2 , q^P^Pa + T 12 (p 2 ,p' 2 , q 2 )p',Jp v p' a 
+T n (p 2 ,p' 2 , q 2 )p^p' u p' a + T u (p 2 ,p' 2 , q^p'^PvPa. (4) 

Due to J? is a conserved current, five constraints among these fourteen independent 
Lorentz structures are introduced. Therefore, only nine of them are independent. However, 
in the sense of calculating coupling constant, we can work with any one of the fourteen 
structures. There are some points that one must follow: (i) The chosen structure must 
appear in the phenomenological side, (ii) The chosen structure should exhibit good OPE 
(operator product expansion) convergence, (iii) The chosen structure should have a sta- 
bility that guarantees a good match between the two sides of the sum rule. After the 
calculations, the structures that obey these points are g^ u p' a in the case p off-shell and 
g^uPa in the case B* off-shell for fB*B* P - Whereas, the structures are g a uP'^ in the case B* 
off-shell, and g afl p' v in the case p off-shell for gB*B* P - 

In order to get the sum rules, the correlation functions need to be calculated in two 
different ways: In phenomenological side, they are presented at the hadron level introducing 
hadronic parameters; in theoretical side, they are calculated in terms of quark and gluon 
degrees of freedom by performing Wilson's OPE. The sum rules for the form factors are 
obtained with both representations being matched via quark-hadron duality and equating 
the coefficient of a sufficient structure from both sides of the same correlation functions. 
In order to improve the matching between the two side of the sum rules, double Borel 
transformation with respect to the variables, P 2 = —p 2 — > M 2 and P' 2 = —p' 2 — > M' 2 , is 
performed. 
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The phenomenological part of the first correlation function is obtained by saturating 
the complete set of appropriate B* and p states. The matrix elements associated with the 
B*B*p momentum dependent vertices can be deduced from Eq. (CQ), which can be written 
in terms of the form factors: 

(B* (p,r ] )p-(q,e)\B*-(p + q,O) = -2V2g B , B , p (q 2 ) (rj* ■ £) (p ■ e*) 

-4V2f B » B » p (q 2 )m B » [( V * ■ e*) ■ O 07* ■ ?)] 

(5) 

The meson decay constants f B * and / p are defined by the following matrix elements: 

(0|jf \B*{p)) = m B *f B *e B *(p), 

m\ P (P)) = m P f P e^(p). (6) 

Saturating Eq. (J2]) with £>* and p states, and using Eqs. (jSJ) and fl6]), then summing over 
polarization vectors via 

^(p)^(p) = -^ + ^, 



e B *{p)e B „*(p) = -g^ + 



m 



2 ' 



(7) 



the physical side of the correlation function for B* off-shell is obtained as 

m p f p f B *m 2 B , 



{P* + m 2 B »){Q* + m 2 B «){P> 2 



mi 



-2VV 



B*B*p\ 



2\/ , ( h> ( lP\( 

rn B * 



a/3 



PaPp,, . PjjC 



m. 



/i f^cB* / 2\ / , ( lv ( lP\( , PaPf\/ . PpPp\ I 



77? 



B* 



7/7 



B* 



+.... (8) 

In a similar way, we obtain the final expression of the physical side of the correlation 
function for an off-shell p meson as: 

m p f p f B *m B « 



(P 2 + m|*)(Q 2 +m2)(P' 2 + 



x 



m 



B* 



2V2g B * B *Jq ){-g^ + -^s )(-#a/3 + —5-) (-^7 + — 



Pa] 



2 ^7 



m B , 7T7.R* m 

„ , PaP-f,, <l, 



-^V2f B *B* P ((l 2 ) m B*(~9 P p + -^)(-#a 7 + ^2 2 )(~^/3 + "2 )li 
1 m B „ m B * mi 



+AV2f B , B , p {q 2 )m B ,{-g„ + ^)(-g a p + 
+..., 



m? B * 



PaPp, ( . QuQp, 



m 2 B 



mi 



(9) 
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B— (e) P B- (f) B'° 

FIG. 1: (a) and (b): Bare loop diagram for the B* and p off-shell, respectively; (c) and (e): 
Diagrams corresponding to quark condensate for the B* off-shell; (d) and (f): Diagrams corre- 
sponding to quark condensate for the p off-shell. 




B— (d) P B— (e) P B— (f) P 

FIG. 2: Diagrams for contributions of bi-gluon operator in the case B* off-shell. 



where represents the contribution of the higher states and continuum. 

In the following, we concentrate our attention to the QCD side of the correlation func- 
tions in the deep Euclidean space. The coefficients Tj above can be written in terms of 
perturbative and condensate terms 

Fi = rf r + rf } + rf } + rf } + r| 6) + • • ■ (10) 

where Tf er is the perturbative contribution, and Ff \ ■ ■ •, F^ are contributions of conden- 
sates of dimension 3, 4, 5, • • • operators in the OPE. The perturbative contribution and 
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FIG. 3: Diagrams for mixed quark-gluon operators in the case B* off-shell, 
gluon condensate contribution can be written in the form of dispersion integration, 



-pper 
p(4) 



4tt 2 
1 



ds 



du 



47T 5 



OO /■OO 

ds I du 



pT>,^,Q 2 ) 

[s — p 2 ){u — p' 2 ) ' 
P?\s,u,Q 2 ) 
s — p 2 ){u — p' 2 ) ' 



14, 



where pi(s,u,Q 2 ) is the spectral density. The spectral density is obtained by calculating 
the bare loop diagrams (a) and (b) in Fig. ([I]) for B* and p off-shell, respectively. In 
the calculation, Cutkosky rules are adopted to deal with the usual Feynman integral 
of these diagrams, i.e., by replacing the quark propagators with Dirac delta function 
2 ^ 2 (— 2iii)5(q 2 — m 2 )6(q°). The integration region for the perturbative contribution 



in Eq. fjllj) is determined from the facts that arguments of the three 5 functions must 
vanish simultaneously. The physical regions of s and u are determined by the following 
inequalities: 



_ i < f b*, s 2su + (s + u-t)(m 2 -s) 

~ l ' ] \^{mls,ml)\^{s,u,t) " ' 
_ (s + u-t)(m 2 + s)-2s(u + m 2 b ) 
~ [ ' >~ \s-m 2 \\V 2 (s,u,t) ~ + ' 



(11) 



where A(a, b, c) = a 2 + b 2 + c 2 — 2ac — 2bc — 2ab and t = q 2 = —Q 2 . 

The quark condensate contribution in the QCD side is determined from the quark 
condensate diagrams (c), (d), (e) and (f) of Fig. ([!]). As what has been shown in Refs. 



16 



191 ]. heavy quark condensate contribution is negligible in comparison with the perturbative 
one. Thus, only light quark condensates contribute to the calculation. It is noticed that 
contributions of diagrams (d), (e) and (f) are zero after the double Borel transformation 
with respect to the both variables P 2 and P' 2 . Hence, we calculate the diagram (c) for the 
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off-shell B* meson and obtain 

n B*(3) m b (uu){-g^ u p' a + g av p'^ + g ail p' v ) 
^ a (p 2 -m 2 b )(p> 2 ) ' 1 } 

The diagrams for the contribution of the gluon condensate in the case B* off-shell are 



depicted in Fig.(T2]). We follow the method employed in Refs. [28|, [29], namely, directly 
calculate the imaginary part of the integrals in terms of the Cutkoskys rule. The diagrams 
for the contribution of the quark-gluon mixing condensate in the case B* off-shell are 
depicted in Fig.Q. The results of the related Borel transformed coefficient BTi in Eq. 
(JTUJ) are given in the appendix. 

The quark-hadron duality assumption is adopted to subtract the contributions of the 
higher states and continuum, i.e., it is assumed that 

p higherstates( s ^ ^ = p OPE^ ^ ^ _ _ ^ ^3) 

where Sq and Uq are the continuum thresholds. 

To improve the matching between the sides of the sum rules, the double Borel trans- 
formation are applied with respect to the P 2 = —p 2 — > M 2 and P' 2 = —p' 2 — > M' 2 . In 
this work we use the following relations between the Borel masses M 2 and M' 2 which are 
= _?|L f or a B * off-shell and = 1 for a p off-shell. 



III. NUMERICAL ANALYSIS 



In the numerical analysis of the sum rules, input parameters are shown in Table [B 
We first determine the three auxiliary parameters, namely the Borel mass parameter M 2 
and the continuum thresholds, sq and uq. The continuum thresholds, sq and uq, are not 
completely arbitrary as they are correlated to the energy of the first excited states with 
the same quantum numbers as the states we concern. They are given by sq = (ms* + A s ) 2 
and uq = (m + A u ) 2 , where m is the p meson mass for the case that B* is off-shell and 
the B* meson mass for that p is off-shell. A n and A s are usually around 0.5 GeV. The 
threshold Sq, Uq and Borel parameter M 2 are varied to find the optimal stability window 
where pole dominance and OPE convergence of the sum rule are satisfied. 



TABLE I: Parameters used in the calculation. 

m 6 (GeV) m B »(GeV) m p (GeV) /s*(GeV) / P (GeV) (mz)(GeV) 3 
4.7 ±0.1 5.325 0.775 0.16 ± 0.01 0.16 ± 0.005 (-0.23) 3 
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(b) 



FIG. 4: a) The OPE convergence of the form factor gB*B*p(Q 2 = 1-0 GeV 2 ) on Borel mass 
parameters M 2 for A s = A u = 0.5 GeV. The notations a, /3, 7, A and p correspond to to- 
tal, perturbative, quark condensate, four-quark condensate and mixed condensate contribution 
respectively and b) pole-continuum contributions. 



A. QCD sum rules for gB*B*p{Q 2 ) 

Using A s = A u = 0.5 GeV for the continuum thresholds and fixing Q 2 = lGeV 2 , 
we show the different contributions to the form factor g B t B * p as a function of the Borel 
variable, as can be seen in Fig. (jlK)- We find an good OPE convergence and a good 
stability of g B *B*p f° r M 2 > 35 GeV 2 . Fig. (jib) demonstrates the contributions from the 
pole term and continuum term with variation of the Borel parameter M 2 . We see that 
the pole contribution is larger than continuum contribution for M 2 < 47 GeV 2 . We choose 
M 2 = 40 GeV 2 reference point. 

Now, we would like to discuss the behavior of the form factors in terms of Q 2 , which is 
plotted in Fig. ([5l). In this figure, the circles correspond to the form factor 9b*b*p{Q 2 ) m the 
interval where the sum rule is valid. Our result is better extrapolated by the mono-polar 
parametrization: 

r>* , ~. 99.4 GeV 2 , N 

= Q* + 75.5 GeV* ' (14) 

Coupling constant is defined as the value of the form factor at Q 2 = —id 



2 where m is 



the mass of the off-shell meson. Using Q 2 = —m 2 B * in Eq. ( 114"]) . the coupling constant is 



obtained as g B *s* P = 2.09. 

In the case p off-shell, Fig. (]6k,) demonstrates a good stability and OPE convergence 
of gs*B*p(Q 2 = 1-0 GeV 2 ) with respect to the variations of Borel mass parameters for 
M 2 > 4 GeV 2 . We see that the pole contribution is bigger than the continuum one in the 
Borel window M 2 < 7 GeV 2 from Fig. ([6b)- We choose M 2 = 6.5 GeV 2 . Our numerical 
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FIG. 5: 5 f 

* B*p(Q 2 ) (circles) QCDSR form factors as a function of Q 2 . The solid line correspond 
to the monopolar parametrization of the QCDSR data. 
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FIG. 6: a) The OPE convergence of the form factor g P B*B* p {Q 2 = 1.0 GeV 2 ) on Borel mass pa- 
rameters M 2 for A s = A u = 0.5 GeV. The notations a, j3 and 7 correspond to total, perturbative 
and four-quark condensate contributions respectively and b) pole-continuum contributions. 



results can be fitted by the exponential parametrization 

-Q 2 



g p B *B* P (Q 2 ) = 0.66 Ex P [ 



0.82 GeV 



2J- 



shown by the solid line in Fig.([7j). Also, g P B*B* P = 1-37 is obtained at Q 2 
(|T5|) . Taking the average of the two results, we get 



(15) 



-m 2 in Eq. 



9b»b- p = 1-73 ±0.25. 



(16) 
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Q [GeV ] 



FIG. 7: g p B 

* B*»( ( 3 2 ) (circles) QCDSR form factors as a function of Q 2 . The solid line correspond 
to the exponential parametrization of the QCDSR data. 



Following the procedure of error estimate in Refs. 30|, |31], with all parameters kept 
fixed, except one which is changed according to its intrinsic error, we calculate a new 
coupling constant and its deviation. Then we obtain percentage deviation related with 
each parameter and how sensitive this value is with respect to each parameter. Table HT1 
show the percentage deviation for the two cases. 





Deviation % 


Parameters 


B* off-shell 


p off-shell 


/b* = 160 ± 10 (MeV) 


14.1 


15.2 


f p = 160 ± 5 (MeV) 


3.2 


7.3 


m b = 4.70 ±0.1 (GeV) 


15.1 


27.8 


M 2 ± 10% (GeV) 


1.8 


1.7 


As ±0.1 e Au±0.1(GeV) 


18.4 


23.5 



TABLE II: Percentage deviation related with each parameter for gs*B' 



p- 



Considering the uncertainties presented in the tables, the coupling constant is: 

9b-b*p = 1-79 ±0.59. (17) 

B. QCD sum rules for fB*B* P (Q 2 ) 



For fB*B*p(Q = 1-0 GeV ), OPE convergence of the sum rule with the Borel mass 
and pole dominance are shown in Fig.®. As the same procedure in the last subsection, 



11 




2 

"a o-i 



C 0.6 

g 
o 

0.4 



20 25 n 30 o 35 40 45 50 20 « 

M 2 [GeVq M 2 [GeV 2 ' 



- Continuum 
-Pole 



(a) 



(b) 



FIG. 8: a) The OPE convergence of the form factor fB*B* P (Q 2 = 1-0 GeV 2 ) on Borel mass pa- 
rameters M 2 for A s = A u = 0.5 GeV. The notations a, j3 and 7 correspond to total, perturbative 
and four-quark condensate contributions respectively and b) pole-continuum contributions. 



the Borel mass is determined to be 25 GeV . With the thresholds A, 



0.5 GeV and 



A u = 0.5 GeV, our numerical calculations of fB*B* P {Q 2 ) can be well fitted by the mono- 



polar parametrization shown in Fig. (Q: 

fE*B* P (Q 2 ) 



13.5 GeV 



Q 2 + 41.6 GeV 2 v ; 

Setting Q 2 = —m 2 Bt in Eq. ( 115]) . the coupling constant is obtained as fs*B* P = 1-01 GeV -1 . 



Figure f|T0|) is plotted to show the OPE convergence, stability and pole dominance 
of the sum rule with the Borel mass for fB*B*p(Q 2 = 1-0 GeV 2 ). With the thresholds 
A s = 0.5 GeV and A u = 0.5 GeV, and a Borel mass of 9 GeV 2 , our numerical calculations 
of fB*B*p(Q 2 ) can be well fitted by the exponential parametrization in Fig. ( ITT]) : 



f B *B* P (Q 2 ) = 0.72 Exp[ 3 44 g e 2 v2 ] GeV -1 , 



(19) 



Using Q 2 = —m 2 in Eq. ( TT9]) . the coupling constant is obtained as fB* B * p — 0.86 GeV . 
Taking the average of the two results, we get 



iWp = (0.94 ± 0.08) GeV" 1 . 



(20) 



Making the same procedure of error estimate as last subsection, Table II III show the 
percentage deviation for the coupling constant b* p - 

Considering the uncertainties presented in the tables, the coupling constant is: 



f B * B * p = (0.94 ± 0.24) GeV" 1 . 



(21) 
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FIG. 9: f§*B*p(Q 2 ) (circles) QCDSR form factors as a function of Q 2 . The solid line correspond 
to the monopolar parametrization of the QCDSR data. 
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(a) (b) 

FIG. 10: a) The OPE convergence of the form factor /b*b*p(Q 2 = 1-0 GeV 2 ) on Borel mass pa- 
rameters M 2 for A s = A u = 0.5 GeV. The notations a, f3 and 7 correspond to total, perturbative 
and four-quark condensate contributions respectively and b) pole-continuum contributions. 



It is noticed that the form factors obtained are different if the B* or the p meson is 
off-shell but both give the compatible coupling constant. As commented in Ref. 16], the 
two sets of points (B* or p off shell) can be fitted by different empirical formulas. However, 
the condition must be satisfied that when extrapolated to Q 2 - 



-m 



where m is the 

mass of the off-shell meson, each fit should go to the compatible value of the coupling 
constant. Together with the predictions from LCSR 13| . the numerical results of the 
coupling constant are presented in Table HT1 Comparison shows that our result gs*B*p and 
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Q 2 [GeV 2 ] 



FIG. 11: fs'B* p(Q 2 ) (circles) QCDSR form factors as a function of Q 2 . The solid line corre- 
spond to the exponential parametrization of the QCDSR data. 





Deviation % 


Parameters 


B* off-shell 


p off-shell 


f B * = 160 ± 10 (MeV) 


15.5 


12.4 


f p = 160 ± 5 (MeV) 


4.5 


3.7 


m b = 4.70 ±0.1 (GeV) 


9.2 


23.4 


M 2 ± 10% (GeV) 


2.8 


1.5 


As ±0.1 e Au±0.1(GeV) 


14.4 


17.9 



TABLE III: Percentage deviation related with each parameter for /b*b* p - 

Jb* B*p are in good agreement with their estimate. However, it is noticed that the authors 
only give the central values of Qb*b* p and fB*B* P , our results are compatible with their 
estimates in case that the uncertainties are considered in their work. 

In summary, the form factors gB*B* P (Q 2 ) and fB*B* P (Q 2 ) parameterizing the B*B*p 
vertex have been calculated in the framework of three-point QCDSR. Both cases that B* 

TABLE IV: Theoretical estimations of the strong coupling constants from different models 

Coupling constant This work [13] 
9B*B* P 1.79 ±0.59 1.88 

/B*B* P ( GeV_1 ) °- 94 ± °- 24 °- 82 
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is off-shell and p is off-shell have been considered. As a side product of the form factors, 
the coupling constants gB*B* P and fB*B* P are estimated, which are compatible with the 
results from the LCSR method 13|. Due to the potential ability in analyzing absorption 
cross sections of T in experiments, the related events are expected to be observed in the 
LHC in the near future. 
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Appendix 

The Appendix is devoted to analytical results of the form factors. Two cases are con- 
sidered, one is for g B * B * p {Q 2 ) and g p B * B * p (Q 2 ), the other for fE* B * P (Q 2 ) and f p B * B * P i.Q 2 )- 
(1) For g$ 

*b* p {Q 2 )i the form factor is 

(QS) = -■">(<? + "30 e ^ e §,BV, (22) 

' /l./ p m|.V2(m|,+mJ + Q') ^ ' 

where 

BY = BY pert + BT {3) + fir (4) + BT {5 \ (23) 
The perturbative contribution is 

BTP ert = -— ds dupf* {pert \s } u } t)6[l-F B *(s,u) 2 }e^e^, (24) 

4 7T Jml JO 

b 

where 

^ rt ^"' i '-pd^)F (s - f+ " ) 

{ml - m 2 b (s + t-u) + st)(u(2m 2 b +s + t)-(s- t) 2 ). (25) 
The condensate contributions are 

Sr (3) = m b (uu)e'^, (26) 

£r (4) = ^- r ds f U ° dupf {4 \ Sj u,t)e[l-F B *(s,u) 2 ]e^e^, (27) 

4 7T Jm? JO 



in which 







2[X(s,u,t)}^ 

((s 3 + 7s(t - uf - lls 2 (t + u) + 3(t - w) 2 (t + m) 
-4m 2 (-2s 2 + s{t - 2u) + (t - u)(t + 2m))), (2* 
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and 



£ r(B) = m b (gqa.G q ) ^ 

AM' 2 v ; 



For g p B * b* P {Q 2 ) ■> t ne form factor is 



where 

BT = BT per + BT (3) + BT {4) + BT {5) . (31) 
The perturbative contribution is 

£ Tper = -l_ r°o ^ ruo duff {p er) {s u t)e[l _ F P (s?M) 2 ]e ^ e ^ ; (32) 
4 7T Z Jml Jml y 



"b "°b 



where 

(ml - m 2 b (s -t + u) + su)(-2m 2 b t + s 2 - s(t + 2u) + u(u - t)). (33) 
The condensate contributions are 

BT {3) = 0, (34) 
BF (4) = f° ds r dup p ^(s,u,t)9[l- F p (s,u) 2 }e^e^ , (35) 

4 IT 2 Jml JO y 

in which 

pf } (s,t,u) 



(g 2 G 2 ) 



2[\(s,u,t)}^ 

(-7s 3 + (t - u) 2 (t + 5m) + s 2 (3t + 7u) + s(3t 2 + 2tu - 5m 2 ) 
+m 2 (8s 2 + 8st + 8t 2 - 4su - 4tu - 4m 2 )), (36) 

and 

£?r (5) = 0. (37) 
(2) For f'B* B * p {Q 2 ), the form factor is 



fZ»J&) = „ ; ~ >Q2 J r f B " ] , . „,A ^BT (38, 



fB*f P m B*rn p ^j2(m 2 B , +m 2 + Q 
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where 

BY = BT per + BT {3) + BT {4) + BT {5) . (39) 
The perturbative contribution is 



[ ds [ dup B ^ per \s,u,t) 

Jm? JO 



A 7T 2 Jm\ JO 

0[1 - F B *(s,u) 2 }e^e 



M' 



(40) 



where 



pf (per \s,t,u) = rw 3U s ,..J 4mlu + mt(-s 2 + 2st - Asu - t 2 - Atu + 5u 2 ) 
1 |A(s, u, t)\ b ' 2 

+m 2 b (2s 2 t + 2s 2 u - Ast 2 + Astu - Asu 2 + 2t 3 - 2t 2 u - 2tu 2 + 2u 3 ) 
-s 2 t 2 - s 2 tu + 2st 3 - 2st 2 u -t 4 + 3t 3 u - 3t 2 u 2 + tu 3 ). (41) 

The condensate contributions are 

BT (3) = 0, (42) 

-1 /•«() fUQ 



Br® = — - ds r dupf (4 \s,u,t)9[l - F B * (s,u) 2 ]e^e^ , (43) 

4 7l 2 Jml Jo 1 



in which 



(9 2 G 2 ) 



2[\(s,u,t)}^ 

(-3s 3 - 3t 3 + s 2 (3t - 7u) - 7t 2 u + lltw 2 - u 3 + s(3t 2 + lAtu + 11m 2 ) 
+4m 2 (s 2 + t 2 + tu + Au 2 + s(-2t + u))), (44) 

and 

BT {5) = 0. (45) 

For fB*B* P {Q 2 )-> ^e form factor is 
where 

BY = BT per + BT {3) + BT {A) + BT {5) . (47) 
The perturbative contribution is 

BTper = ^ dupf per \s,u,t)6[l - F p (s,u) 2 }e^e^, (48) 

4 7T Z Jml Jml J 
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where 

pf per \s,t,u) = * (t(-2mg(a + t-u)+ mg(4s 2 - 2a(t + n) - 2(t - m) 2 ) 

-m 2 (s + t - u)((s - tf + 4sw - 2tw + ii 2 ) + 2su(s(t + u)-(t- u) 2 ))). 

(49) 

The condensate contributions are 

5r (3) = 0, (50) 
BT^ = ^\ r ds r dup p f (4 \s,u,t)9[l- F p (s,u) 2 }e^eTfi , (51) 

4 Tl 2 Jm\ Jm\ 1 

in which 



pP/(4)(M ' M) " 2[A(s,M)]5/2 

(-9s 3 + {t- u) 2 (5t -u) + s 2 (llt + 5m) 

+M fe 2 (8s 2 + 8st + 8t 2 - 4sm - 4tw - 4m 2 ) + s(-7t 2 + lOtu + 5m 2 )), (52) 

and 

5r (5) = 0. (53) 



^ 2 G 2 ) 
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